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q ' Abstract 

The effect of the conducting electrode on the interaction of dust 
^ (— i grains in a an ion flow is discussed. It is shown that two grains levitat- 

Qh ing above the electrode at the same height may attract one another. 

This results in the instability of a dust layer in a plasma sheath. 

> _ 

^ ■ 1 Introduction 

o 

A good deal of experiments in dusty plasma physics are performed with 
■ aerosol grains suspended in the plasma sheath. A negatively charged dust 

grain in the sheath levitates above a horizontal electrode due to the balance 
of two forces: the gravity force directed downwards and the sheath electric 
• ^ \ field that pushes the grain upwards. The near-sonic or the supersonic ion flow 

in the sheath creates the wake field downstream of the grain. Since the latter 
is confined by a certain Mach cone, it is commonly accepted that the wake 
field affects the motion of grains, which are situated downstream only. The 
usual assumption is that the intergrain potential is smooth in the horizontal 
direction, i.e., two grains levitating at the same height repel one another via 
the screened Coulomb potential. The structure of the wake field and the 
grain interaction in an ion flow were studied in details [1-5]. The asymmetric 
interaction of vertically aligned grains was also observed experimentally ||. 
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Analytical theory and computer simulations cited above assumed that 
the plasma density is constant and the influence of the conducting electrode 
is negligible. Although both assumptions evidently fail under conditions of a 
real plasma sheath, taking into account the plasma non-uniformity seems to 
be a very difficult problem. In order to estimate the influence of the electrode 
upon grain interaction, here I use the zeroth approximation, which however, 
seems the only one treatable analytically. 

In the present paper, the following simple model is accepted. Let there 
be the monoenergetic ion stream entering a conducting electrode (or a grid) 
located at the horizontal plane, z — 0. The stream velocity, u, exceeds the 
ion thermal velocity but it is much less than the electron thermal velocity. 
Electrons are Boltzmann distributed. Two problems are addressed: first, how 
does the electrode modify the interaction between two grains levitating at 
the same height and, second, how this affects the spectrum of dust acoustic 
waves propagating along a single dust layer. 

2 Intergrain interaction 

The electrostatic potential produced by a point charge, Q, located at r = r 
is given by the solution of the Poisson equation 

Ae>(r) = -47rQ<y(r-ro), (1) 

where e is the operator of the static dielectric permittivity of an ambient 
plasma. Within the accepted model, the spatial Fourier transform of e is 
given by 

eft) = 1- 71 XW + f> (2) 
(k z u + i0) 2 k l 

where u is the velocity of the ion flow, Ui is the ion plasma frequency and 
ko is the inverse electron Debye length. The ion flow is parallel to the z axis 
and directed downwards, u < 0. 

In unbounded plasma, the natural boundary condition for Eq. (|TJ) is 
<f\r-^oo — ► 0. It is convenient to express the solution to Eq. flU) in terms 
of its Green function, i.e., ip(r) = QG°(r — tq), where the Fourier transform 



2 



of G°(r) with respect to the transverse coordinate = (x,y) is 

Here kj_ = (k x , k y ) stands for the transverse components of a wave vector. 
With the dielectric function given by Eq. (0), the zeros of denominator in 
Eq. (D are 

k z = ±k D q — iO, 

k z = ±ik D n, (4) 

where 




q± = k±/kjj is the normalized transverse wave vector, and /i = Ui/kr>u is the 
inverse Mach number. 

The integral in Eq. (§) is readily evaluated resulting in 

G IM = TT-r—2 [*e~"* ,M +29(-z)q sin (qk D z)] . (7) 
ko q + k 

The spatial structure of the potential is recovered by the Fourier transform 
of the Green function (0) with respect to k±. The first term in parenthesis in 
expression (|7|) gives rise to the Debye-Hiickel potential distorted by the ion 
flow, while the second term represents the wake field situated downstream of 
the charge [1-3]. The opening of the Mach cone confining the wake and the 
field structure inside it depend on the stream velocity. 

Now we turn to the evaluation of the electric potential of a charge located 
near a conducting wall. Let the wall be situated at z = plane, while the 
charge is placed above it, Zq > 0. Then the potential is given by the solution 
of Eq.(|l|) supplemented with the boundary condition (p\ z= a = 0. As it is 
well-known from electrostatics, we can make allowance for this boundary 
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condition by introducing the surface charge density, aj™^ induced at the 
conducting surface. Then, the potential of a unit charge is written as 

G kx (z, z ) = G° k± (z - z ) + G° k± (z)al n ± d . (8) 

Taking into account the boundary condition, G k± (0,z ) = 0, we find the 
surface charge density a l k nd and, finally, 

n i \ - no ( \ G l ± ( z ^ G k±(~ z o) / n \ 
^k±{z,z ) — Lr k± {z - z ) G°~(0) ' 

One may doubt whether description a bounded dispersive medium in 
terms of the response function of an unbounded medium is justifiable. How- 
ever, more accurate and lengthy calculations give the same result. The phys- 
ical reason is that there are no ions reflected by the conducting wall within 
the present model. The mathematical reason is that Eq. ([]]) actually masks 
a set of partial differential equations with two real characteristics directed 
downwards. 

Of particular interest for the following is the interaction potential of two 
charges placed at the same height, z$. Since 

4tt e ~ Ka 

G k± (zq, z q ) = — — - [k sinh(rea) + q sin(ga)] , (10) 

kd k + q 

where a = koZo, the normalized potential is 

oo 

/ e — Ka 
q±dq± Jo(q±p) 2 2 sinh(Ka) + gsin(ga)] , (11) 
K + q 

o 

where G(r±, zq, Zq) = kdiv(p) and p = kdr. The asymptotic behaviour of 
w(p) is determined by Eq. ( P7G|) at qj_ — > 0. The latter depends essentially 
on whether the the ion flow is supersonic or subsonic. In the case of the 
supersonic flow, p < 1, the roots (|]||) at q± — » are approximated as 

q^q±^==, K W yfl -p 2 , (12) 
V 1 -M 2 

and the leading term of the asymptotic expansion of w(p) (O) is 



P \/a l p l - p 2 (l - p 2 ) 
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This expression may be interpreted as a mirror reflection of the wake field 
produced by the grain at < z < z . More detailed numerical investigation 
of the potential (|T3| ) shows that w(p) is always positive for \i < 1. 

Quite another behaviour is observed for the case of the subsonic flow, 
fi > 1. The roots (|5|||) at q± — > are now 

q « vV - 1, K « Q± 7== , ( 14 ) 
while the potential at infinity behaves like 

w(p)\ p ^oo ~ 2°^ 1 Sin (qV)" 2 ~ l) ( 15 ) 

The most important distinction between the expressions (|T3|) and (|T5|) is 
that in the subsonic regime the potential is attractive if 



sin a 



vV - i) < 0. (16) 



The numerically evaluated example of the potential fll5|) demonstrating 
the attraction is depicted in Fig. [|. It should be pointed out that the in- 
equality flTBD guarantees the long-scale attraction between grains. With the 
opposite inequality imposed on a and //, the potential is repulsive at large 
distance but the attractive branch may appear at smaller scales, as shown in 
Fig-©- 

Evidently, the existence of the attractive branch of the interparticle in- 
teraction may result in formation of various patterns and clusters even in 
the absence of the external potential well confining grains in the horizontal 
direction. Also, the even distribution of grains in a dust layer may become 
unstable. The latter possibility is discussed in the next section. 



3 Dust layer 

Now consider a two-dimensional gas consisting of dust grains hovering over 
the conducting electrode. Ignoring intergrain correlations, the linearized 
equations of motion are written as 
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Figure 1: The distribution of the potential in the transverse direction, fi = 2, 
a = 2.7. 




Figure 2: The distribution of the potential in the transverse direction, /i — 2, 
a = 4. 
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JI + V ± (a v) = 0, (17) 
^ = -uv-^k d V ± Jdv' ± w(kD(v ± -v' x )), (18) 

where <t is the unperturbed value of the surface density and a is the density 
perturbation. Here I consider horizontal motions only, i.e., v z = 0. The term 
— w in Eq. flTS] ) corresponds to the grain friction on an ambient neutral gas. 
The intergrain interaction is described by the potential given by Eqs. (p"T|). 
Although the grain charge, Q, generally depends on the ambient plasma 
parameters, for simplicity we ignore its variability. 

Assuming that all quantities are proportional to exp(— ioot + ik±r±) we 
easily get the dispersion relation for the gas oscillations: 

k 2 

uj(uj + iv) = g d - L G k± (z ,z ), (19) 
lix 

where g d = 2irQ 2 a /M. 

This expression describes dust sound waves in the continuous medium 
approximation. In the long- wave limit, q± — > 0, the dispersion relation is of 
the form 



+ iu)= f dql \*-*^^*V^), M<1 (20 ) 
V |1 — M 2 | [sm(a v //i 2 - 1), /i>1 

Evidently, the layer is unstable, i.e., Imu > 0, if Gk ± (z , z ) < 0. In the long- 
wave limit this is possible in the subsonic flow only and the corresponding 
constraint coincides with Eq. flip]). More detailed investigation shows that the 
potential, ^^(^0,-20)5 is always positive if \i < 1. However, in the subsonic 
regime, \i > 1, there are regions of instability, that is, Gk ± (z ,z ) < if 
< q± < qmax(a, /•*)■ The latter are shown as shadowed areas in Fig. ^. With 
the increasing distance (a — > 00) to the wall or decreasing stream velocity 
(/i — > 00) the instability regions shrink to zero, q m ax(a, ^) ~ a -1 / 2 , fi^ 1 ^ 2 . 



4 Conclusion 

To summarize, we have shown that the presence of the conducting wall may 
drastically change the electrostatic interaction of the dust grains in an ion 
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Figure 3: Instability regions in the (q±,a) plane, /z = 2. 

flow. In particular, the electrostatic image of the grain wake field may result 
in attraction between grains levitating at the same height that, in its turn, 
yields the Jeans-type instability of the dust layer. 

It would be unduly naive to draw quantitative conclusions from the 
present calculations. However, it seems reasonable that even in a real plasma 
sheath, which is essentially non-uniform, the electrostatic image of the grain 
wake field may also affect the motion of another grains outside the Mach 
cone. Although in most experiments the screened Coulomb interaction is 
observed, there are indications that the intergrain potential may be more 
complicated. It was recently reported that, under certain conditions, a void, 
i.e., a dust-free region, appears in a central part of a single dust layer J7J. The 
emergence of a two-dimensional void in a layer consisting of some hundreds 
of grains can hardly be explained in the same manner as a three-dimensional 
void; the latter requires strong influence of the dust component upon the 
discharge structure [§]. Although currently one cannot exclude that some 
additional external forces appeared in the experiment 0, we can conjecture 
that the 2D void formation is provided by complicated intergrain interaction, 
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for example, the one described in this paper. 
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